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Optimal  Replacement  Age  in  an 
Imperfect  Inspection  Model 

by 

Donna  C.  Herge,  Frank  Proschan 
and  Jayaram  Sethuraman 


ABSTRACT 

A  device  is  maintained  under  an  age  replacement  policy.  The  status  of  the 
device  (functioning  or  failed)  is  known  only  by  inspection  at  some  fixed  interval 
k.  With  probability  q,  an  inspection  error  may  be  made,  and  a  functioning  unit 
will  be  declared  to  have  failed  and  be  replaced  by  a  new  unit.  On  the  contrary, 
when  a  failed  unit  is  inspected,  it  is  assumed  that  no  inspection  error  will  be 
made.  Assuming  that  the  cost  of  replacing  a  failed  unit  (actually  failed  or 
believed  failed)  is  greater  than  the  cost  of  replacing  a  functioning  unit,  we 
show  that  we  can  obtain  an  optimum  replacement  age  which  minimizes  L(T),  the 
long-run  expected  cost  per  unit  of  time.  We  find  a  lower  bound  for  the  optimal 
replacement  age  and  obtain  asymptotic  and  monotonicity  properties  for  L(T). 
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1 .  INTRODUCTION 


Most  of  the  recent  work  on  maintenance  models  takes  into  account  additional 
real-world  factors  like  "minimal"  or  "imperfect”  repair  actions.  Another  relevant 
factor,  imperfect  inspection,  has  not  been  considered  extensively  in  the  literature. 

Brown  and  Proschan  (1983)  propose  several  different  imperfect  maintenance 
and  imperfect  inspection  models  in  [3]  and  study  properties  of  the  distribution 
of  the  mean  time  between  perfect  repairs  in  [4],  Fontenot  and  Proschan  (1984) 
then  develop  optimal  policies  for  several  maintenance  models  based  on  the  imperfect 
repair  model  in  [4], 

The  objective  of  this  paper  is  to  develop  an  optimal  replacement  policy  for 
an  imperfect  inspection  model.  Derman  and  Sacks  (1960)  solved  the  problem  of 
choosing  an  optimal  replacement  rule  for  deteriorating  equipment  when  the  amount 
of  deterioration  is  observed  periodically.  Perfect  inspection  was  assumed. 

Barlow,  Hunter,  and  Proschan  (1963)  showed  how  to  obtain  optimum  inspection 
schedules  for  a  broad  class  of  failure  distributions  in  [1].  For  an  unknown 
failure  distribution,  Derman  (1961)  obtained  the  minimax  schedule  assuming  that 
failure  has  probability  p>0  of  being  detected  by  an  inspection,  and  in  the 
absence  of  failure  the  unit  is  replaced  at  some  specified  time  T.  More  recently, 
Taylor  (1975),  assuming  a  cumulative  damage  model  for  system  failure,  found  an 
optimal  replacement  strategy  that  minimizes  the  long-run  expected  cost  per  unit 
of  time. 

In  this  paper  we  are  concerned  with  a  device  subject  to  an  age  replacement 
policy.  We  assume  that  all  repair  (replacement)  actions  are  perfect  and  repair 
time  is  negligible.  The  state  (functioning  or  failed)  of  the  device  is  determined 
by  periodic  inspection  at  a  specified  interval  k.  The  failure  of  a  unit  remains 
undetected  until  the  unit  is  actually  inspected.  Due  to  an  inspection  error 
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occurring  with  probability  q,  a  functioning  unit  is  declared  to  have  failed  are 
is  replaced  by  a  new  unit.  On  the  contrary,  when  a  failed  unit  is  inspected,  no 
error  in  inspection  occurs. 

An  inspection  error  may  be  due  to  human  error  or  malfunction  of  a  detection 
device.  Inspection  policies  are  frequently  used  for  safety  or  security  devices. 
One  example  is  a  fire  extinguisher  which  is  checked  periodically.  If  the  pressure 
gauge  malfunctions,  it  may  erroneously  indicate  the  extinguisher  is  empty  and 
so  the  extinguisher  is  unnecessarily  replaced. 

It  is  apparent  that  an  age  replacement  policy  is  inappropriate  if  the  under¬ 
lying  failure  rate  is  decreasing.  Thus  we  will  assume  that  the  failure  rate 
r(t)  =  f(t)[F(t)]-1  is  increasing.  We  assume  that  the  life  distribution  F  of 
the  device  is  absolutely  continuous  with  density  f  and  F(0)  =0. 

We  use  the  notation  F  =  1  -  F,  F^  =  F(i),  and  r^  =  r(i). 

2.  FORMULATION  AND  SOLUTION  OF  MODEL. 

A  device  has  life  length  X  with  distribution  F.  The  device  is  installed 
at  time  0,  and  is  inspected  at  successive  times  k,  2k,  ...  .  When  the  device  is 
functioning,  the  probability  of  no  error  in  inspection  is  p  =  1  -  q.  The  device 
is  replaced  by  a  new  unit  at  age  T,  which  is  a  positive  integer  multiple  of  k, 
or  at  the  first  inspection  following  failure,  whichever  comes  first.  This  process 
is  continued  indefinitely.  Clearly  the  process  renews  itself  at  times  of 
replacement. 

The  probability  that  the  device  is  replaced  after  time  t  =  ik  is  P[Y  >ik]  * 
p1  for  i  =  0,  1,  2,  ...,  since  the  device  must  survive  to  time  ik  and  pass 
i  inspections.  Thus  Y  is  the  observed  life  length  of  the  device.  Note  that 
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EY  -  ^  I  F.,  p1  .  Let  Z  be  the  elapsed  time  between  replacements;  then 
i=0  1K 

Z  =  min(Y,T). 

Let  Cj  be  the  cost  of  an  unscheduled  replacement  (actual  or  believed 
failure)  and  be  the  cost  of  a  scheduled  (age)  replacement.  Assume  c^  <  C2  . 

Our  objective  is  to  find  the  replacement  age  T  which  minimizes  L(T) ,  the  long-run 
expected  cost  per  unit  of  time.  In  Theorem  2.1  we  give  a  formula  for  L(T)  and 
in  Theorem  2.2  we  describe  the  optimal  value  of  T.  These  results  are  analogous 
to  those  found  by  Fontenot  and  Pros chan  (1984)  for  a  modified  age  replacement 
model  with  imperfect  repair. 


Theorem  2.1 


L(T) 


( cl"c2^ ET-k 


Proof:  The  times  of  replacement  are  renewal  times  in  3  renewal -reward  process 
whose  interarrival  distribution  is  that  of  Z  =  min(Y,T).  From  renewal  theory. 


we  have  that 


LCD 


.  C(T) 

“  dHT 


> 


where  C (T)  is  the  expected  cost  per  renewal  cycle  and  P(T)  is  the  expected 
duration  of  a  renewal  cycle.  Denote  the  distribution  function  of  Y  by  G.  Then 


5(t)  =  P[Y>t]  =  for  (i-l)k  s  t  <  ik 

and  i  =  1,  2,  ... 
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Then  the  expected  length  of  a  cycle  is 


T  T/k  ik 

D(T)  =  E(Z)  =  /  G(t) dt  =  l  /  6(t)dt 

0  i=l  (i-l)k 


T/k  .  ,  ik  IT-1  . 

-  I  ^fi-nk  p  /  dt  =  kI  PikP  • 

i=l  U  1Jk  (i-l)k  i=  0  1K 


With  probability  F_,  k  p  ,  the  device  is  replaced  after  time  T  -  k  and  a  cost  of 

c2  is  incurred.  Otherwise,  it  is  replaced  at  or  before  time  T-k  and  a  cost  of 

Cj  is  incurred.  Therefore  the  expected  cost  per  cycle  is 
T  T 

r--l  ^-1 

cm  =  c1(l-FT_kpK  )  +  c2  FT_kP  .  o 


Theorem  2.2.  In  addition  to  the  assumptions  about  F  in  Section  1,  assume  that 
r  is  differentiable.  Then 

T 

T_i  £-1 

a.  LTH  =  F_  ,pF  [k  l  F..  pi]'2H(T),  where 


H(T)  =  (cj-c^  [(krT_R- logp)  ?ik  p1  +  FT,k  pX  ]  -  Cj. 

b.  An  optimal  replacement  age  T*  exists  (T*  may  be  infinite)  and  T*  =  k  if  and 

c2 

only  if  p  s  exp(- — —  +  krn) . 

c2-<-i  u 

cl  1 

c.  T*  is  finite  if  and  only  if  rro  >  fc'  ygy  +  Flogp* 

£Y  C1 

d.  —  <  —  implies  that  T*  <  “  . 
k  c2 

c2 

e.  A  lower  bound  for  T*  is  — -  EY  . 

C1 


4 


▼ 


Proof:  a.  Differentiate. 

b.  Given  >  c7  note  that  H  is  increasing  since  r  is  increasing  and 

H'(T)  =  (cj-cp  k  r'_k  jo  F.^1  . 

Then  T*  =  k  <=>  L  is  increasing  <=>  H(k)  s  0  <=>  (e^-^Mkrg-  logp  +  1)  _^c^  £  0 
C2 

<=>  p  <  exp(- — — -  +  kr  ). 

c2-cl  u 

c2 

If  p  >  exp( - +  krn),  then  H(k)  <  0.  Thus  if  we  can  find 

c2“cl  u 

Tq  =  min{T:  H(T)SO),  then  L“  is  negative  on  (0,Tg-k)  and  positive  on  (TQ>®). 
So  we  choose  T*  such  that  L(T*)  =  min{L(T0-k) ,L(Tq) }.  Otherwise,  if  there  is 
no  T  such  that  H(T)  s  0,  then  L  is  decreasing  and  T*  =  <*>. 

c.  The  inequality  involving  r  implies  lim  H(T)  >  0.  Thus  the  equation 

T-x® 

TQ  =  min{T:  H(T)>0}  has  a  finite  solution. 

c2  C1  C2 

d.  Note  that  L(T*)  S  L(k)  =  /k.  If  T*  =  ”,  then  L(T*)  =  /EY  <;  4/k.  By 

EY  C1 

the  contrapositive,  /k  <  / Cj  implies  that  T*  <  °°. 

°2  C1 

e.  It  is  easy  to  show  that  L(T)  2:  .  Thus  (LT)  must  cross  L(°°)  =  —  (if  it 

1  c  c  EY 

''l 

crosses  at  all)  to  the  right  of  the  value  at  which  —  crosses  —  ;  this  value 

1  fcY 

‘■'2 

is  T  =  — EY.  Therefore,  we  should  schedule  replacement  at  an  age  greater 

than  —  EY.  □ 

C1 

3.  Asymptotic  and  monotonicity  properties  of  L(T). 

If  we  let  pel,  we  get  a  desired  asymptotic  result  for  L(T)  as  the  length  of 
the  inspection  interval,  k,  goes  to  zero.  In  this  case  our  model  reduces  to  the 
simple  age  replacement  model  for  which  the  optimization  problem  was  solved  by 
Barlow  and  Proschan  (1965),  pp.  85  -  90.  Since  the  denominator  of  L(T)  is  a 
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Pieman  sum,  it  converges  to  the  integral  form  given  in  Barlow  and  Proschan  (1965) 


p.  88,  if  F  is  Riemann  integrable. 


lim  L(T)  =  lim 
k-*0  k-K) 


^2  +  ^cl“c2^  pT-k  c2  +  (ej_c2)F,p 

k-K)  X-l  T  - 

k  f  F(x)dx 

7  k  F.,.  0 


l 


If  0  <  p  <  1  we  will  show  that  L(T)  as  k  -►  o.  This  result  is  expected 
since  more  frequent  inspections  will  lead  to  more  units  being  erroneously  replaced. 


Case  1.  If  T  =  k,  then  lim  L(k)  =  lim  -r-  =  «. 

k->0  k-K)  K 


ci  ci  ci 

Case  2.  If  T  =  «,  then  lim  L(«)  =  lim  =  lim - = -  >  lin - - — 

k"°  k*°  W  kl  F^p1  k^°  k  l  p1 

i=0  i=0 


*  lim 


*  *  ***  «  « 

k-K)  kf-U 
'■l-p-' 


Case  3.  If  k  <  T  <  °°,  then  lim 

k-»-0 


cl-(crc2)FT.kP 
k  I  PiiP1 


>  lim 


C1  ‘  ^cl_c2^  p 


k-0  ~-l 

k  I  p1 

i=0 


=  lim 


C1 *  (ci-c2>  P 


k<0  k(^-) 

-  b 1 

For  a  particular  value  px  of  p,  we  denote  L  (T)  =  [c1  -  1 

T  -  Pi 

F_1 

^ikpi  *  Then  it:  is  easily  seen  that  Lp(T)  increases  as  p  f  0  as  follows. 
i=0 


a.  Note  that  F^p*  s  F-k  p2x  for  Pj  <  p2,  i  -  0,  1,  ...,  £  -  1. 
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'..11011  T>k,  we  have  stric‘ 


T  T 

--1  ir-1 


IT*  .  k 

Thus  1/ [k  l  l/[k  l  fiyP^  ]  for  T  =  k,  2k,  ...  . 


i='. 


i  =  0 


ineaualit.y  above. 

T  T 

£-1  _  1-1 

b.  Note  that  FT  k  Pj  £  p*  for  Pj  <  p2  and  T  =  k,  2k,  ...  .  When  T  >  k 


b> 


h1 


we  have  strict  inequality.  Thus  c^  -  (°i"c2^x-k  P1  ^  C1  ~  ^cl”c2^T-k  p2 

From  a  and  b  we  have  Lp^(T)  =  Lp^(T)  for  T=k  and  Lp^ CT)  >  Lp^(T)  for 
Pj  <  p2  and  T  =  2k,  3k,  ...  .  This  is  intuitively  reasonable  since  more  mistakes 
in  inspection  lead  to  an  increase  in  unnecessary  replacements.  This  can  be  seen 

in  Figure  4.10  in  the  next  section.  Clearly  lim  L  (T)  =  c./k  . 

pkO  P 

Since  Lp(T)  is  a  decreasing  function  of  p,  we  obtain  a  sharp  lower  bound  for 
L(T)  by  taking  the  limit  as  p  increases  to  1.  Thus 


'-i  “  (''i“c2^T-k 

L(T)  ^  lim  L  (T)  =  — -  X-  — -  =  L  (T)  . 

pH  P  F-l 

k  I  P-V 
i=0  lk 

This  bound  is  shown  in  Figure  4.10. 

We  denote  the  optimal  value  of  T  for  a  particular  p  as  T*.  Then  the  optimal 

long-run  expected  cost  per  unit  of  time  is  L  (T*)  =  min  L (T) .  Since  for  all  T 

P  P  t  p 

we  have  Lp^(T)  £  Lp^(T)  ^or  P1  <  p2’  **  ^°^ows  t^iat  =  myn  2 

min  l_,  (T)  =  L_  (T*  ).  Thus  L  (T*)  is  also  a  decreasing  function  of  p.  This 
T  p2  p2  p2  P  P 

can  be  seen  in  figure  4.10. 


4.  APPLICATION’S 

Our  model  makes  sense  only  for  a  failure  rate  which  is  increasing.  We  see 
from  Theorem  2.1  that  the  denominator  of  L(T)  is  a  sum  which  can  be  explicitly 
solved  for  an  exponential  distribution  F.  Although  it  is  clear  that  the  optimal 
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policy  is  to  replace  only  at  failure  times,  because  replacing  an  exponential 
component  by  a  new  component  does  not  change  the  residual  life,  we  will  show 
that  Theorem  2.2  confirms  that  T*  =  «  .  Let  F(x)  =  e_Ax,  X  >  0,  and  k  =  1. 

oo 

Then  EY  =  £  (^t)1  =  — - -  .  Setting  L'(T)  =  0  and  solving  for  T  we  get 

i=0eA  l-P/eX 


T 


1  +  log 


2 


1  + 


)  ( 1--^- )  +  log  \ 

e  e 


Note  that  1  +  (log  -E--1)  >  o  since  it  reduces  to  the  form  logx<x  -  1,  which 

6  X  e  c 

is  true  for  x  =  >  1.  Thus  for  T  to  exist,  we  need  ( — — )  (l-iL)  +  iog-£-  >  0, 

P  C.-C-,  ex  gA 

which  implies  - — —  >  —■".,*??  P  f  which  in  turn  implies  r  =  X  <  7 - ~Ycv  +  logp. 

l“c2  l-p/eA  "  Cc1-c2JbY 

Thus  from  Theorem  2.2c.  we  see  that  T*  =  ®. 

Next  we  will  consider  an  application  of  our  model  to  the  truncated  normal 

distribution  so  that  we  can  compare  our  results  to  those  in  Example  1,  p.  90,  of 

Barlow  and  Proschan  (1965).  "Example  1.  Many  types  of  electron  tubes  used  in 

commercial  airline  communication  equipment  and  elsewhere  tend  to  have  a  truncated 

normal  failure  distribution  (Aeronautical  Radio,  Inc.,  1958).  A  certain  tube 

used  in  commercial  equipment  has  a  truncated  normal  failure  distribution  with  a 

mean  life  of  9080  hours  and  a  standard  deviation  a  of  3027.  Suppose  Cj^  =  $1100 

and  c2  =  $100,  so  that  c^/c2  =  11." 

The  density  f(x)  of  the  truncated  normal  distribution  may  be  written  as 


f(x) 


0  otherwise. 


A 
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where  *(x)  =  — —  e'x  /2  and  b  *  -  7  *  ( 2^)  dx  .  Fro*  L'(T)  =  0  we  get 

a  0  0 

T 


-1 


(2.1) 


Ckr^-logp)  l  Plk  P1  .  FT.kpr 


-1 


i*0 


crc2 


Making  the  change  of  variable  y ^  and  defining  r^(x)  =  —  ,  we  have 


from  (2.1)  for  a  truncated  nomal  distribution  F, 

oy0>ii 


fl>(  t)dt 
x 


oy0+v 


(2.2)  K(y0)  *  TN(y0)  -  logp)  Ip1/  *(v)dv  +  p  k  /  <Hv)dv 

i=0  ik-y  yQ 

a 

be,  be- 

1  b  +  £ 


C1"C2 


cl“c2 


If  2:  3,  then  f(x)  is  very  close  to  the  density  of  a  normal  distribution  with 
mean  u  and  standard  deviation  o. 

In  Example  1,  Barlow  and  Proschan  estimate  b  by  1  and  use  a  graph  of 

K(yQ)  =  K(yg)  +  b  for  ^  *  3  to  find  y^.  To  obtain  more  accurate  results  for 

bc_ 

comparison  with  our  model,  we  compute  b  =  .9987  and  K(y.)  »  - - a  .09987.  Then 

rc2 

Yq  =  -1.63  and  the  optimal  replacement  age  is  4146  hours  with  an  associated 
minimum  cost  of  1,(4146)  =  $.036. 

For  our  model  we  consider  several  examples  with  various  values  for  k,  the 

length  of  the  inspection  interval,  and  p,  the  the  probability  of  no  error  in 

inspecting  a  functioning  unit.  Note  that  for  the  truncated  normal  distribution, 

r®  =  °°»  so  w®  know  by  Theorem  2.2c.  that  T*  is  finite.  Using  (2.2)  we 

T-k-u  be, 

wish  to  find  the  value  of  T  such  that  K(—  —  )  is  nearest  to  K(yn )  =  - - 

O  ' Q  C\~C2 

1.0986. 


i 

1 
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Using  (2.2)  we  compute  the  following  values: 


T 

y 

K(y) 

1000 

-3.00 

1.050 

2000 

-2.67 

1.051 

3000 

-2.34 

1.060 

4000 

-2.01 

1.088 

5000 

-1.68 

1.149 

FIGURE  4.1:  Values  of  K(y)  for  k  =  1000,  p  =  .95  . 


We  see  from  Figure  4.1  that  T*  =  4000  is  the  optimal  value,  which  gives  us 

the  minimum  cost  L(T*)  =  $.0710  as  shown  in  Figure  4.2.  We  compute  EY  =  7629. 

Recall  from  the  proof  of  Theorem  2.2e  that  c^/EY  =  .144  is  the  horizontal 

c2 

asymptote  for  L(T).  Also  —  EY  =  694  is  a  lower  bound  for  T*. 

c! 

We  see  from  Figure  4.3  that  T*  =  4000  is  the  optimal  value  and  L(T*)  =  $.1271 

c2 

shown  in  Figure  4.4.  We  compute  EY  =  5979.  Then  —  EY  =  544  is  a  lower  bound 

c,  c.  L1 

for  T*.  Notice  in  Figure  3.4  that  >  gy  *  t*ius  implication  in  Theorem  2. 2d 
cannot  be  reversed. 

Comparing  Examples  2  and  3,  it  is  clear  that  more  frequent  inspections 
result  in  a  higher  optimal  cost  L(T*). 

We  see  from  Figure  4.5  that  the  graph  of  K(y)  does  not  cross  K(yQ)  =  1.0986. 

Thus  there  is  no  solution  to  L^T)  =  0,  so  that  L'  >  0,  and  the  cost  is  strictly 

increasing,  as  shown  in  Figure  4.6.  Then  the  minimum  cost  is  $.10  =  L(1000)  at 

T*  =  k.  This  is  the  expected  result  from  Theorem  2.2b  since 
c~ 

.5  =  p  S  exp( -  +  krn)  *  exp[-.l  +  1 000 (.00000 15 ) ]  «  .9062.  Since  EY  =  1980, 

C2“cl  U 

c2 

— EY  =  180  is  a  lower  bound  for  T*. 

C1 

Just  as  in  Example  4  we  see  from  Figure  4.7  that  K  does  not  cross 

K(y^)  *  1.0986,  so  that  cost  is  strictly  increasing  and  T*  =  k  as  expected  from 

c_ 

Theorem  2.2b  since  .5  =  p  £  exp  ( -  +  krn)  =  exp [-. 1  +  500 (.0000015) ]  x  .9055. 

C2'C1  U 

Q 

A  lower  bound  for  T*  is  ~EY  =  (998)  =  91. 

Cj  1100 


4500 


Example  4: 


1000 


T 

y 

K(y) 

1000 

-3.00 

1.6923 

2000 

-2.67 

1.5410 

3000 

-2.34 

1.4710 

4000 

-2.01 

1.4490 

5000 

-1.68 

1.4590 

6000 

-1.35 

1.4969 

FIGURE  4.5:  Values  of  K(y)  for  k  =  1000,  p  =  .5  . 
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.50 
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HOURS 


4000 
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Example  5:  k  =  500,  p  =  .5 


T 

y 

KCy) 

500 

-3.00 

1.6915 

1000 

-2.83 

1.5384 

1500 

-2,67 

1.4630 

2000 

-2.50 

1.4260 

2500 

-2,34 

1.4096 

3000 

-2.17 

1.4040 

3500 

-2,01 

1.4042 

4000 

^1“ 

00 

0 

1 

1.4089 

4500 

-1,68 

_ 

1.4164 

FIGURE  4.7:  Values  of  K(y)  for  k  =  500,  p  =  .5. 


EXAMPLE  5 


L(T) 


1 

EY 


1 . 1  0  2_ 
1.10  % 


1.00 


.90 

.80 

.70 

.60 

.50 

.40 

.30 


o 

o 

m 


o 

o 

o 

ro 


o 

O 

m 

m 


o 

o 

o 

J- 


*  T 


HOURS 


4500 


We  summarize  the  results  as  follows: 


EXAMPLE 

P 

k 

c2 

—  EY 

C1  P 

tp 

EY 

P 

L  (T*) 

P  P 

1 

1.00 

0 

c2 

~U=825 

C1 

4146 

p=9080 

$.036 

2 

.95 

1000 

694 

4000 

7629 

$.071 

3 

.95 

500 

544 

4000 

5979 

$.127 

4 

.50 

1000 

180 

1000 

1980 

$.100 

5 

.50 

500 

91 

500 

998 

$.200 

Table  4.9:  Summary  of  Examples. 

Asymptotic  and  monotonicity  properties  of  L^CT)  are  shown  in  Figure  4.10 
for  k=  1000.  Optimal  values  Lp(T*)  are  circled.  For  p=  .5,  .7,  and  .9  we  have 
T*  «  k  *  1000.  Clearly  L^CT*)  and  L^CT)  are  decreasing  in  p.  Lj(T)  is  a  sharp 
lower  bound.  Horizontal  asymptotes  are  labeled  Lp(°°).  Lines  connecting  values  of 
Lp(T)  for  a  given  value  of  p  were  drawn  to  illustrate  the  trend  of  Lp(T)  as 
T  -*■  ®. 
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REPLACEMENT  AGE,  T  HOURS 


Asymptotic  and  monotonicity  properties  of  Lp(T)  for  an  inspection 
interval  k  =  1000  hours,  for  a  truncated  normal  life  distribution. 
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